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It is shown that a given non-autonomous system of two first-order ordinary differential equa- 
tions can be expressed in Hamiltonian form. The derivation presented here allow us to obtain 
previously known results such as the infinite number of Hamiltonians in the autonomous case 
and the Helmholtz condition for the existence of a Lagrangian. 
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Se muestra que un sistema dado, no autonomo, de ecuaciones diferenciales ordinarias de primer 
orden puede expresarse en forma hamiltoniana. La deduccion presentada aquf nos permite 
obtener resultados previamente conocidos tales como el mimero infinito de hamiltonianas en el 
caso autonomo y la condicion de Helmholtz para la existencia de una lagrangiana. 
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1. Introduction 



As is well known, it is very convenient to express a given system of ordinary differential equations 
(not necessarily related to classical mechanics) as the Euler-Lagrange equations associated with 
some Lagrangian, L, or as the Hamilton equations associated with some Hamiltonian, H (see, 
e.g., Ref. 1). One of the advantages of such identifications is the possibility of finding constants 
of motion, which arc related to symmetries of L or H. Also, the Hamiltonian of a classical 
system is usually regarded as an essential element to find a quantum version of the mechanical 
system. 

In the simple case of a mechanical system with forces derivable from a potential (that may 
depend on the velocities), there is a straightforward procedure to find a Lagrangian or a Hamil- 
tonian. However, in the case of non-conservative mechanical systems or of systems not related 
to mechanics, the problem of finding a Lagrangian or a Hamiltonian is more involved. A given 
system of n second-order ordinary differential equations are equivalent to the Euler-Lagrange 
equations for some Lagrangian if and only if a set of conditions (known as the Helmholtz con- 
ditions) are fulfilled (see, e.g., Refs. 2, 3 and the references cited therein). 

The aim of this paper is to give a straightforward procedure to find a Hamiltonian for a 
given system of two first-order ordinary differential equations (which may not be equivalent to a 
second-order ordinary differential equation) that possibly involves the time in an explicit form. 
The results derived here contain the Helmholtz condition for n = 1 (in the case where the given 
system is equivalent to a second-order equation). In Sec. 2 the main results of this paper are 
established, demonstrating that a given system of first-order ordinary differential equations can 
be expressed in Hamiltonian form looking for an integrating factor of a differential form made 
out of the functions contained in the system and, in Sec. 3, several examples are presented. In 
Sec. 4 we show that, in the appropriate case, our results lead to the Helmholtz condition for the 
existence of a Lagrangian. 

2. Hamiltonians and canonical variables 

We shall consider a system of first-order ordinary differential equations of the form 

x = f{x,y,t), y = g{x,y,t), (1) 
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where / and g are two given functions. A system of this class can be obtained from a second-order 
equation 

X = F{x, X, t), 

by introducing the variable y = x. We are initially interested in finding a Hamiltonian function, 
H, and canonical variables, q, p, such that the corresponding Hamilton's equations be equivalent 
to the system ([T|). 

Assuming that there is an invertible relation between the variables x, y and a set of canonical 
coordinates q, p, x = x{q,p,t), y = y{q,p,t), in such a way that Eqs. ([T]) are equivalent to the 
Hamilton equations for q and p with a Hamiltonian H, making use of the chain rule, one finds 
that 

d(x^ yj dy dx 

— gdx + fdy = —-^-—dH — ;:;-dx + —dy + terms proportional to dt. (2) 
o{q,p) ot at 

Therefore, given the system ([1]) we start by considering the differential form 

-{g- 0)dx + (/ - V)dy, (3) 

where 

dy{q,p,t) dx{q,p,t) 
nq,P,t) = g-^ , '^[Q,P,t) = — 

are functions unspecified by now (see Eq. ()10p below). For a fixed value of t, the differential 
form ([3]) is always integrable (see any standard text on ordinary differential equations, e.g., Ref. 
4); that is, there exist (locally) functions a and H, which may depend parametrically on t, such 
that 

-{g- 4>)dx + (/ - i,)dy = adH. (4) 

Now, for simplicity, without any loss of generality (since, once we have found a set of canonical 
coordinates, we have the liberty of making any canonical transformation afterwards), we choose 
q = x (hence, ^ = 0) and, therefore, 

d{x,y) dp 
d{q,p) dy' 

Then, by comparing Eqs. ^ and the canonical momentum, p, must be such that 

dp{x,y,t) 1 



dy a{x,y,t)' 

Hence 



(5) 



dp 1 dp 
dp = —dx + -dy + —dt (6) 
dx a dt 



or, equivalently, 



dp dp 
dy = —a—dx + adp — a—dt (7) 



thus, recalling that x = q, this last expression shows that 

dp{x,y,t) 
'^ = -^^^ 

and we can also write Eq. ([6]) in the form 



(8) 



dv 1 d> 
dp = -i-dx + -dy - -dt. (9) 
ox a a 

Since this is an exact differential, we have 

This equation establishes a relation between the integrating factor and the function (see 

examples below). 

From Eqs. dH), with -i/' = 0, and ([9|) we have 

1 1 f)}-f 

dH = —{g- 0)dx + -fdy + —dt 
a a at 



■-(g- ^)Ax + / f dp - -^Ax + j + —it 

»_* + /|P)d, + /dp+('§ + /*)d*. 

a a dx J \ dt a J 



Hence, considering H function of q, p, and t, 



[see Eqs. ([T])] and 



dH g (p dp . 

= + f^=P, (12) 

dq a a dx 



since, according to Eqs. ([9]) and ([T]), 

dp y (j) dp g cj) 

p= + = — /H . 

dx a a dx a a 

Equations (llip and p2p are equivalent to the original system ([1]) and have the desired Hamil- 
tonian form. 

Summarizing, the system of equations ([1]) can be written in the form of the Hamilton equa- 
tions, with the Hamiltonian determined by Eq. ^ and the canonical momentum defined by Eq. 



The fact that the left-hand side of Eq. ^ multiphed by a ^ is an exact differential yields 
(when ip = 0) 

|;[--fa-*)l = ^(.-'/), 



which amounts to 



b-0)^ + .-'-f(i,-« + /^ + .-'|^=O. (13) 

ay ay ox ox 



Hence, making use of Eqs. ([T]), (fT3|) and ([TO]) , we obtain 



d _i da da da 
-T7<^ = — — X H — y + 



dt dx dy dt 

da~^ da~^ da~^ 



(?f + ¥) ■ (14) 



dx dy dt 

. da-' _ ^-1 _ 0) _ ^-1^ + ^ 
dy dy dx dt 

'dl^dg_^ 

^dx dy y 

(Note the cancelation of (/>.) 

Equation (|14|) shows that the function a is determined up to a factor that is a constant of 
motion and, therefore, there exists an infinite number of Hamiltonians (and, correspondingly, of 
expressions for p). It may be noticed that Eq. ([M]) is just Liouville's theorem. 

3. Examples 

A first example is provided by the equation 

X + 7X + oJqX = r]{t), 

where 7 and loq are constants, and rj{t) is an arbitrary function, which corresponds to a forced 
damped harmonic oscillator. Taking y = x, we have y = —^y — oJqX + r/(t), which is of the form 
([1]) with /(x, y, t) = y, and g{x, y, t) = —^y — uJqX + r]{t). Then Eq. (fT^ reduces to 

d -1 -1 
d^" 

and we can take a = e"'*'* (any other choice would require the knowledge of the explicit form of 
r]) then from Eq. pO|) we see that 

del) 
dy 



which is satisfied with (j) = —^y. Substituting all these expressions into Eq. (jH) we have (with t 
treated as a constant) 



and, therefore, we can take H = e^\y^ /2 + uolx^ /2 - r]{t)x). Finally, from Eq. dlD we find that 
p can be chosen as p = e^^y. The corresponding Lagrangian can be calculated in the usual way, 
by means of the Legendre transformation. 

The results of the previous section allow us to readily derive those of Ref. 5, corresponding to 
the autonomous case. In fact, when the functions / and g, appearing in Eqs. (1), do not depend 
explicitly on the time, from Eqs. ^ and ([T]), taking (p = = ip, we have aH = —gx + fy = 
—gf + fg = 0. This means that H is some constant of motion, which is not unique; we can 
replace it by H' = G{H), with G being an arbitrary function. H' is also a constant of motion and 
a will not depend explicitly on t [see Eq. (llOp ]. no matter what (time-independent) Hamiltonian 
we choose. 

The expressions given above allow us to find H, which need not be related to the total energy. 
In the example considered in the appendix of Ref. 5, f{x,y) = y, g{x,y) = —ky, where /c is a 
constant (i.e., x = —kx). Then, —gdx + fdy = kydx + ydy = yd{kx + y) and, therefore, we can 
take CT = y and H = kx + y. 

We end this section by considering the problem studied in Ref. 6 (which corresponds ap- 
proximately to a relativistic particle subjected to a constant force. A, and a force of friction 
proportional to the square of the velocity) , namely (with the appropriate changes in notation) 



where m represents a mass. A, 7, and a are constants. Thus, f{x,y) = y, and g{x,y) = 
{X--fy^){l-a^y^)/m. Thus, 



{ujqX — r]{t))dx + ydy = e ^^dH 



my = {X- 7/)(l - a^y^) 



-gdx + fdy = (A - 72/^) (1 - a^y^)dx + ydy 




Comparing with Eq. @ (with (p = = ip) we immediately see that we can take 



a = {X - ^y^){l - aV) 



and 




6 



X 1 

m ~^ 2(Aa2 — 7) 



1 — a^y^ 

According to Eq. the canonical momentum p can be taken as 

dy 



(A-7y2)(l-aV)' 

Despite the huge difference with the expressions given in Ref. 6, one can show that the Hamil- 
tonian obtained in that reference is essentially the exponential of our H. (See, Eqs. (23) and 
(26) of Ref. 6.) 

4. The Helmholtz condition 

The case where one starts with a second-order equation of the form 

x = F{x,x,t) (15) 

(considered in Refs. 2, 3), is a particular case of the treatment above if one defines, e.g., y = x, 
that transforms Eq. (|15p into the system 

x = y, y = F{x,y,t), 

which is of the form ([T|) with f{x, y,t) = y and g{x, y, t) = F{x, y, t). Then Eq. reduces to 

which is the Helmholtz condition when there is one degree of freedom (see, e.g., Ref. 2 and the 
references cited therein; note that cr~^ = dp/dy = dp/dx = d'^L/dxP' is the integrating factor 
wii employed in these references). 

On the other hand, not every system of equations of the form ([1]) comes from a second-order 
equation x = F{x, x, t). An example is given by 

x = f{x,t), y = g{y,t). 



where there is no coupling between the variables x, y. Here (choosing (/> = = V') 

dx dy 

T^7 



-gdx + fdy = fg ( + 
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Therefore, if we assume that cr = fg does not depend explicitly on of t [see Eq. pO|) ]. we can 
take 



Thus, H = pf — J f ^dx and with the Hamiltonian being a linear function of p, the Legendre 
transformation is not defined nor the Lagrangian. 

5. Concluding remarks 

As we have shown, at least in the case of a system of two first-order ordinary differential equa- 
tions, finding a Hamiltonian is essentially equivalent to finding an integrating factor for a linear 
differential form in two variables. The integrating factor also determines the expression for the 
canonical momentum. Equation (jl4p is analogous to the Helmholtz condition, but, in the present 
approach, it leads directly to the Hamiltonian (in the standard approach, finding a solution to 
the Helmholtz conditions, only gives the second partial derivatives d'^L/dxidxj). When the 
system is non-autonomous, it is convenient to find the integrating factor using Eq. (|14p . while 
in the autonomous case, it may be more simply obtained from the linear differential form itself. 
Finally, as shown in Sec. 4, there are systems of equations for which a Lagrangian does not exist, 
but a Hamiltonian description can be given. 
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